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Abstract — WindSat  is  the  world’s  first  microwave  polarimetric 
radiometer,  designed  to  measure  ocean  vector  winds.  In  late 
2004,  the  first  preliminary  oceanic  wind  vector  results  were 
released,  and  this  paper  presents  the  first  evaluation  of  this 
product  for  several  Atlantic  hurricanes  during  the  2003  season. 
Both  wind  speed  and  wind  direction  comparisons  will  be  made 
with  surface  wind  analysis  (H*Wind)  developed  by  the  NOAA 
Hurricane  Research  Division  (HRD)  and  provided  by  the 
NOAA  National  Hurricane  Center  (NHC).  Examples  are 
presented  where  HRD  aircraft  flights  were  conducted  within 
several  hours  of  the  WindSat  overpass.  These  H*Wind  surface 
wind  analyses  provide  the  most  complete  independent  surface 
winds  comparison  data  set  available.  Both  WindSat  retrieved 
wind  speeds  and  wind  directions  are  evaluated  (against 
H*Wind)  as  a  function  of  storm  quadrant.  To  complement  the 
analysis,  rain  rates  were  derived  using  WindSat  brightness 
temperatures  with  a  modified  version  of  the  TMI  2A12 
heritage  rain  algorithm.  Effects  of  rain  on  the  derived  wind 
speeds  and  directions  are  discussed. 

I.  Introduction 

WindSat  is  the  first  space-bome  microwave  polarimetric 
radiometer,  launched  on  the  Coriolis  near-polar  LEO  satellite 
in  January  2003.  As  a  proof-of-concept  mission,  the 
objective  of  WindSat  is  to  develop  and  demonstrate  passive 
ocean  wind  vector  retrieval  from  space  as  risk  reduction  for 
the  upcoming  Conical  Scanning  Microwave  Imager/Sounder 
(CMIS),  planned  for  launch  in  late  2009.  Following  the 
heritage  of  seven  Special  Sensor  Microwave/Imagers 
(SSM/I),  WindSat  provides  retrievals  of  atmospheric  and 
ocean  environmental  parameters  such  as  total  cloud  liquid 
water,  integrated  water  vapor,  wind  speed  and,  additionally, 
sea  surface  temperature.  The  focus  of  WindSat,  however,  is 
to  provide  the  first  passive  wind  direction  measurements 
from  space.  Currently,  two  preliminary  six-month  data  sets 
of  WindSat  retrievals  are  available:  an  empirically  derived 
product  from  NOAA/NESDIS  [1,2]  and  a  physics-based 
product  from  the  Naval  Research  Lab  [3].  Both  products 


demonstrate  that  passive  wind  vector  retrievals  are  feasible 

[4]. 

Ocean  vector  wind  data  is  vital  in  understanding  global 
climate  change,  air-sea  interaction  and  atmospheric 
circulation.  Of  special  interest  to  meteorologists  is  the  ability 
of  space-borne  instruments  to  sense  tropical  cyclones  (TC)  in 
openocean  where  other  instruments,  such  as  buoys,  are 
unavailable.  While  scatterometer  measurements,  including 
those  from  QuikSCAT,  have  been  introduced  into  numerical 
models  used  for  forcasting  hurricanes,  there  has  only  been 
limited  success  due  to  the  detrimental  effects  of  high  rain  on 
scatterometer  retrievals. 

The  objective  of  this  paper  is  to  analyze  the  suitability  of 
passive  ocean  wind  vector  measurements  in  TCs.  Version-0 
wind  retrievals  produced  by  National  Oceanographic  and 
Atmospheric  Administration  National  Environmental 
Satellite,  Data  and  Information  Service  (NOAA/NESDIS) 
are  compared  to  an  experimental  surface  wind  analysis, 
H*Wind  for  two  passes  from  2003.  The  correlation  of  wind 
retrieval  error  with  rain  is  also  investigated. 

II.  Passive  Microwave  Polarimetry 

The  microwave  ocean  brightness  temperature  (Tb)  scene 
observed  by  a  radiometer  depends  upon  various  atmospheric 
and  ocean  surface  geophysical  variables,  as  well  as  a  number 
of  instrument  parameters.  Since  the  Tb  signatures  are 
functions  of  the  radiometer  frequency  (wavelength)  and 
polarization,  the  separation  of  these  effects  allow  for  multi- 
frequency  radiometer  geophysical  retrievals  process  [5].  Of 
special  interest  are  the  signatures  of  ocean  wind  vectors. 
Recent  advances  in  polarimetric  radiometry  modeling  and 
measurements  have  demonstrated  that  usable  wind  direction 
information  can  be  obtained  by  combining  the  vertical  and 
horizontal  polarizations  with  the  cross-correlation  of  the  two 
[6].  The  cross  correlation  terms  represent  the  third  and 
fourth  parameters  of  the  modified  Stokes  vector  [7],  defined 
as 
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In  this  definition,  Tv,  Th,  T+45,  T_45,  Tic  and  Trc  represent 
brightness  temperature  (radiances)  at  vertical,  horizontal, 
plus  45°,  minus  45°,  left-hand  circular  and  right-hand  circular 
polarizations  respectively.  In  principle,  the  Stokes  vector 
provides  a  full  characterization  of  the  electromagnetic 
signature  of  the  ocean  surface,  which  is  sufficient  to  uniquely 
determine  the  wind  direction. 

The  3rd  (implemented  by  the  difference  between  the  ±45° 
polarization  channels)  and  the  4th  (implemented  by  the 
difference  between  the  left  and  right-hand  circular 
polarization  channels)  Stokes  parameters  offer  the  most  wind 
direction  information.  On  WindS at,  there  are  three  channels 
there  are  three  polarimetric  channels  available:  10.7  GHz, 
18.7  GHz  and  37  GHz.  Both  the  3rd  and  4th  Stokes  signatures 
are  relatively  weak  as  environmental  radiometric  signals  go; 
but,  fortunately,  they  are  very  robust  in  that  they  are  highly 
immune  to  the  influences  of  atmospheric  absorption  and 
emission  caused  by  water  vapor  and  cloud  liquid  water. 
These  atmospheric  signals  are  common-mode  and  cancel 
during  the  subtraction  of  the  ±45°  and  left  and  right-hand 
circular  polarizations.  Unfortunately  the  effect  of 
heterogeneous  precipitation  is  less  certain  because  there  are 
an  extremely  limited  number  of  observations  of  the  3rd  and 
4th  Stokes  parameters  in  the  presence  of  rain.  Since  strong 
bands  of  precipitation  accompany  hurricane  winds,  an 
analysis  of  the  effects  of  rain  on  hurricane  retrievals  is 
presented  in  this  paper. 

III.  Data  Sets 

A.  WindSat  Wind  Vector  Retrievals 

This  paper  contains  an  analysis  of  the  WindSat 
NOAA/NESDIS  version-0  wind  vector  retrieval  algorithm 
[1,2]  within  TCs.  Since  Tb  channels  are  of  different 
resolution  and  geolocation,  brightness  temperatures  are 
resampled  to  every  fourth  37  GHz  location  and  beam 
averaged  to  the  resolution  of  6.8  GHz  channel,  resulting  in  a 
spatial  resolution  of  approximately  50  km  at  12.5  km  spacing 
in  both  the  along-  and  cross-track  directions.  The  first  step  in 
this  wind  retrieval  process  is  to  process  the  four-frequency, 
vertically  and  horizontally  polarized  WindSat  Tb’s  in  a 
similar  manner  to  the  procedure  implemented  for  SSM/I  [8] 
to  solve  for  four  relevant  environmental  variables  (not 
including  wind  direction).  A  forward  radiative  transfer 
model  provides  the  relationship  between  the  measured  Tb’s 
(over  oceans)  and  a  variety  of  environmental  parameters.  Of 
these  environmental  parameters,  the  four  that  can  be 
retrieved  are  sea  surface  temperature  (SST),  surface  wind 
speed  (WS),  integrated  atmospheric  water  vapor  (WV),  and 
integrated  cloud  liquid  water  (CLW).  Using  regression 
coefficients  derived  from  the  forward  model,  the  retrieval 
algorithm  sequentially  solves  first  for  the  CLW  (using  the 


18.7,  23.8  and  37.0  GHz),  then  WV  (using  same  channels), 
next  WS  (10.7,  18.7,  23.8  and  37.0  GHz),  and  finally  SST 
(using  only  10.7  GHz).  After  resolving  these  four 
parameters,  the  entire  process  is  iterated  several  times  until 
the  retrievals  converge  to  a  stable  value.  It  should  be  noted 
that  this  environmental  retrieval  algorithm  does  not  include 
the  effects  of  rain,  nor  has  it  been  developed  for  winds  >  20 
m/s. 

The  retrieved  wind  speeds  are  then  used,  along  with  the 

10.7,  18.7  and  37  GHz  3rd  and  4th  Stokes  parameters,  for 
wind  direction  retrieval.  Relative  wind  directions  are  varied 
and  modeled  3rd  and  4th  Stokes  parameters,  (TU5y)mod  are 
calculated.  For  each  relative  wind  direction  trial,  the  square 
of  the  difference  between  the  measured  and  the  modeled 
parameters  {(TU;V)meas,  and  (TU;y)mod}  is  calculated  and 
summed  for  all  frequencies.  Values  of  relative  wind  direction 
that  minimize  the  sum  of  these  terms  for  all  frequencies  are 
“possible  solutions”.  Due  to  the  harmonic  nature  of  the  3rd 
and  4th  Stokes  parameters,  retrievals  results  in  multiple 
possible  wind  directions  (called  wind  direction  aliases  or 
ambiguities).  Aliases  are  ranked  inversely  proportional  to 
the  cost  function.  After  ranking,  the  aliases  subjected  to  a 
multi-pass  median  filter  to  produce  a  single  selected  vector. 
The  alias  selection  skill,  which  is  defined  as  the  percentage 
of  the  selected  directions  that  are  closest  to  the  true  direction, 
is  typically  >  85%. 

B.  WindSat  Derived  Rain  Rates 

To  aid  in  the  evaluation  of  hurricane  wind  retrievals,  we 
derive  rain  rates  using  a  preliminary  WindSat  rain  algorithm. 
The  algorithm  is  based  on  a  modified  version  of  the  Goddard 
Profiling  Algorithm  (GPROF)  currently  being  used  with  the 
Advanced  Microwave  Scanning  Radiometer  -  Earth 
Observing  System  (AMSR-E)  [8,9].  GPROF  is  a  multi¬ 
channel,  physically  based  rainfall  rate  retrieval  algorithm 
employing  cloud-resolving  models  to  produce  a  large 
number  of  cloud  profiles  for  various  rain  conditions. 
Brightness  temperatures  are  produced  using  a  forward  model 
and  are  then  convolved  with  the  antenna  patterns  of  the 
various  frequencies  generating  modeled  Tbs  for  each  of  the 
utilized  channels.  Once  the  database  of  profiles  and 
respective  Tbs  is  established,  rainfall  rates  are  estimated,  on 
a  pixel-by-pixel  basis,  using  a  Baysian  inversion,  given  the 
set  of  measured  brightness  temperatures. 

For  rain  retrieval,  WindSat  brightness  temperatures  are 
resampled  to  every  fourth  37  GHz  location  using  an  average 
of  the  four  nearest  neighbors  with  a  1/R2  weighting.  This 
effectively  yields  a  12.5  km  by  12.5  km  grid,  on  which  each 
channel  retains  its  inherent  resolution.  Four  frequencies  and 
seven  channels  are  used:  10.7,  18.7,  37  GHz  (vertical  and 
horizontal),  and  23.8  GHz  (vertical  only).  Typically  GPROF 
uses  both  37  GHz  and  high-resolution  85  GHz  brightness 
temperatures  to  classify  rain  events  as  convective  or 
stratiform.  Since  the  85  GHz  channel  is  not  available,  the 
WindSat  algorithm  depends  only  on  the  emission 
characteristics  of  the  37  GHz  channel  to  assess  the  type  of 
rain  event  [10]. 


C.  H*Wind  Modeled  Wind  Fields 

Due  to  the  poor  temporal  and  spatial  sampling  of 
hurricanes,  verification  of  TC  surface  winds  is  inherently 
difficult;  however,  with  data  available  from  airplane 
reconnaissance  missions,  reasonable  analyses  are  possible. 
The  H*Wind  [11]  analysis  system  is  a  software  tool, 
developed  by  NOAA’s  Hurricane  Research  Division  (HRD), 
that  assimilates  these  aircraft  measurements,  along  with 
surface  measurements  and  remotely  sensed  satellite  data  into 
a  common  model  that  allows  for  limited  human  quality 
control.  There  are  some  caveats  associated  with  using 
H*Wind  data.  H*Wind  is  not  an  instantaneous  view  of  a 
TC,  but  an  assimilation  of  data  over  a  period  of  three  to  six 
hours.  Also,  much  of  the  storm  it  not  observed  over  this  time 
frame.  Due  to  the  high  temporal  and  spatial  variability  of  TC 
wind  field  structures,  there  is  a  significant  error  range — on 
average  10%  to  20%  for  wind  speed  [12].  Additionally,  the 
under-sampling  of  TCs  results  in  maximum  wind  estimates 
that  are  typically  lower  than  the  actual  intensity. 

IV.  Analysis 

For  this  paper,  two  cases  from  2003  are  carefully  chosen 
where  the  TCs  have  been  well  sampled  by  air  reconnaissance 
close  to  the  respective  WindSat  passes:  Fabian  on  3 
September  2147  UTC  and  Isabel  on  17  September  1040 
UTC.  Using  National  Hurricane  Center  (NHC)  “best-track” 
data,  the  centers  of  the  H*Wind  analyses  have  been  shifted  in 
time  to  match  the  storm  centers  at  WindSat  pass  time.  To 
match  the  spatial  resolution  and  geolocation  of  WindSat 
wind  vectors  and  H*Wind  fields,  a  nearest-neighbor 
interpolation  is  performed,  averaging  all  H*Wind  data  within 


a  25  km  of  radius  of  each  WindSat  measurement.  Since  the 
objective  of  this  paper  is  to  determine  the  potential  of 
WindSat  to  make  reasonable  wind  direction  estimates 
hurricane,  and  not  to  evaluate  selection  skill,  all  analysis 
compares  the  WindSat  directional  alias  that  is  closest  to  the 
H*Wind  direction. 

Fig.  1  shows  a  visual  comparison  between  the  H*Wind 
wind  field  (a)  and  the  WindSat  retrieval  (b)  for  Hurricane 
Fabian  3  September  2147  UTC.  WindSat  wind  speeds  in  the 
northern  region  of  the  storm  are  obviously  biased  high,  while 
retrievals  to  the  south  of  the  eye  are  biased  low.  If  we  look 
at  Fig.  4,  we  see  that  the  overestimates  of  wind  speed 
correlate  well  with  areas  of  intense  rain.  Again,  from  an 
objective  analysis,  wind  direction  retrievals  also  seem  to 
suffer  from  rain  contamination.  While  WindSat  does  seem 
to  represent  the  overall  circulation  of  the  TC  well,  there  are 
some  regions,  even  without  rain  contamination,  where 
WindSat  does  a  poor  job  of  resolving  the  general  flow,  such 
as  on  the  southeast  side  of  the  storm.  In  Fig.  5,  a  quantitative 
analysis  demonstrates  a  high,  almost  linear,  correlation  of 
error  wind  speed  (a)  and  direction  (b)  versus  rain  rate.  RMS 
errors  for  speed  and  direction  are  6.64  m/s  and  29.7° 
respectively. 

V.  Conclusion 

Using  a  two  cases  from  the  2003  hurricane  season,  this 
paper  takes  a  first  look  at  evaluating  the  performance  of 
WindSat  hurricane  wind  vector  retrievals.  Comparisons  of 
wind  speeds  demonstrate  that  the  version-0  algorithm  does 
not  perform  well  for  high  wind  and  high  rain  conditions. 
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Figure  1.  Hurricane  Fabian  2003  3  September  2147  UTC.  (a)  is  H*Wind  surface  analysis,  (b)  is  WindSat  wind  retrieval. 
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Figure  3.  Hurricane  Fabian  2003  3  September  2147  UTC.  WindSat 
wind  directions  superimposed  over  WindSat  rain  rates. 

While  the  performance  of  wind  directions  is  slightly  better, 
there  still  are  major  deficiencies  in  WindSat’s  ability  to 
measure  wind  direction  in  the  presence  of  rain.  We  are 
encouraged,  however,  that  WindSat  is  able  to  make  retrievals 
in  heavy  rain.  Previous  radiometers,  such  as  SSM/I  do  not 
make  wind  speed  retrievals  in  the  presence  of  rain,  and 
traditional  scatterometer  wind  vector  retrievals  in  hurricane 
conditions  are  highly  suspect.  Since  the  version-0  algorithm 
was  developed  only  for  rain  free  wind  retrievals  at  <  20  m/s, 
we  feel  that  further  algorithm  development  will  only  improve 
WindSat’ s  ability  to  retrieve  wind  vectors  in  turbulent 
weather  conditions,  such  as  those  found  in  hurricanes. 
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Figure  5.  Scatter  plots  of  wind  speed  error  (a)  and  wind  direction  error  (b)  as  a  function  of  rain  rate.  X  represents  mean  error  for  each  2.5  mm/hr  bin  of  rain 

rate.  Error  bars  denote  +/-  1  standard  deviation.  Diagonal  is  linear  regression. 


